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ABSTRACT 

We study the clustering properties of the recently compiled SDSS cluster cat- 
alog using the two point correlation function in redshift space. We divide the 
total SDSS sample into two richness subsamples, roughly corresponding to 
Abell R > and APM clusters, respectively. If the two point correlations 
are modeled as a power law, £(r) = (r /V) 7 , then the best-fitting parame- 
ters for the two subsamples are r Q — 20.7+3'g hr 1 Mpc with 7 = 1.6ig '\ and 
r = 9.7+}'2 with 7 = 2.0+Q'g ft." 1 Mpc, respectively. Our results are consistent 
with the dependence of cluster richness to the cluster correlation length. 

Finally, comparing the SDSS cluster correlation function with that pre- 
dictions from three flat cosmological models (Q m = 0.3) with dark energy 
(quintessence), we estimate the cluster redshift space distortion parameter 
P ~ tt^/bo and the cluster bias at the present time. For the ACDM case we 
find P = 0.2±g^g 
scale cluster motions. 

Keywords: galaxies: clusters: general - cosmology: theory - large-scale 
structure of universe 
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1 INTRODUCTION 



Galaxy clusters occupy a special position in the hier- 
archy of cosmic structure formation, being the largest 
gravitationally collapsed objects in the universe. There- 
fore, they appear to be ideal tools for testing theories of 
structure formation as well as studying large-scale struc- 
ture. The traditional indicator of clustering, the cluster 
two-point correlation function, is a fundamental statis- 
tical test for the study of the cluster distribution and 
is relatively straightforward to measure from observa- 
tional data. 

Indeed, many authors based on optical and X-ray 
data have shown that the large scale clustering pat- 
tern of galaxy clusters is well described by a power law, 
£ (r) = (r /r) 7 , with 7 = 1.6 — 2. The correlation length 
r lies in the interval r = 13 — 25/i _1 Mpc, depending on 
the cluster richness as well as the analyzed sample (cf. 
Bahcall & Soneira 1983; Klypin & Kopylov 1983; La- 
hav et al. 1989; Bahcall & West 1992; Peacock & West 
1992; Dalton et al. 1994; Nichol, Briel & Henry 1994; 
Croft et al. 1997; Abadi, Lambas & Muriel 1998; Bor- 
gani, Plionis & Kolokotronis 1999; Collins et al. 2000; 
Tago et al. 2002; Moscardini, Matarrese & Mo 2001; 
Gonzalez, Zaritsky & Wechsler 2002). However, a seri- 
ous issue here is how the galaxy clusters trace the under- 
lying mass distribution. The cluster distribution traces 
scales that have not yet undergone the non-linear phase 



of gravitationally clustering and thus simplifying their 
connection to the initial conditions of cosmic structure 
formation. Galaxy clusters is strong biased with respect 
to the matter distribution (e.g. Peacock & Dodds 1994 
and references therein). 



In this paper we utilize the recently completed 
SDSS CE cluster catalog (Goto et al. 2002) in order: 
(i) to study the two point correlation function in red- 
shift space and (ii) to calculate the relative cluster bias 
at the present time comparing the observational results 
with those derived from three flat cosmological models 
with dark energy (quintessence). The structure of the 
paper is as follows. The observed dataset and its mea- 
sured correlation function are presented in section 2. In 
section 3 we give a brief account of the method used to 
estimate the predicted correlation function in different 
CDM spatially flat cosmologies. The linear growth rate 
of clustering in quintessence cosmological models can be 
found in section 4, while in section 5 we fit the SDSS 
cluster clustering to different cosmological and biasing 
models. Finally, we draw our conclusions in section 6. 
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Figure 1. The estimated (histogram) and the expected (line) 
number of the SDSS clusters as a function of redshift. 



2 ESTIMATION OF THE SDSS CLUSTER 
CORRELATION FUNCTION 

2.1 Cluster catalogue 

In this work we use the recent SDSS CE cluster cata- 
log (Goto et al. 2002), which contains 2770 and 1868 
galaxy clusters in the North (145.1° < a < 236.0°, 
-1.25° < 8 < 1.25°) and South (350.5° < a < 56.61°, 
— 1.25° < S < 1.25°) slices respectively, covering an 
area of ~ 400deg 2 in the sky. Redshifts are converted to 
proper distances using a spatially fiat cosmology with 
Ho = 100/ikms" 1 Mpc" 1 and fi m = 0.3. The cluster 
redshifts are estimated using the color information by 
identifying the bin in g — r which has the largest num- 
ber of galaxies around the color prediction of elliptical 
galaxies (Fukugita, Shimasaku, Ichikawa 1995) at dif- 
ferent redshifts (which define the different g — r bins). 
Due to the fact that the true and estimated redshifts 
are better correlated for z < 0.3 (Goto et al. 2002), we 
will limit our analysis within this redshift range, corre- 
sponding to a limiting distance of r max < 836/i -1 Mpc. 

In Fig. 1, we present the estimated (histogram) and 
the expected for a volume limited sample (solid line), 
number of the SDSS clusters as a function of redshift. 
It is evident that the number of SDSS clusters appears 
to follow the equal-volume oc r 3 law out to z ~ 0.23, a 
fact corroborated also from the standard Kolmogorov- 
Smirnov (KS) test which gives probability of consistency 
between model and observations (up to z < 0.23) of 
Vks — 0.43. Therefore, this SDSS cluster sample is the 
only to-date sample that is volume limited to such a 
large distance and can thus play an important role in 
large scale structure studies. 

We apply the cluster correlation function analy- 
sis using clusters of two richness class: (a) N sa ,\ > 
30 members (roughly corresponding to Abell R > 
0; hereafter Si sample) and (b) iVg a i > 20 mem- 
bers (roughly corresponding to APM clusters; here- 
after S2 sample). These two subsamples contains 200 
and 524 entries with corresponding mean densities of 
n Sl (< 2 m „) ~ 8.42(±0.06) x Hn 6 h 3 Mpc _3 and ns 2 (< 
z max ) ~ 2.20(±0.10) x 10 _5 ft 3 Mpc~ 3 , giving rise to 
inter-cluster separations of the order of ds 1 ~ 49.15 ± 
2.56 /i _1 Mpc and ds 2 ~ 35.66 ± 2.18 /i _1 Mpc, respec- 
tively. 
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Figure 2. The spatial two-point correlation function 
(points) in redshift space for the Si (Abell R=0 richness- 
right panel) and S2 (APM richness-left panel) samples. The 
error bars are estimated using the bootstrap procedure. The 
dashed lines represent the best-fitting power low £gj ( r ) 
'?' (see parameters in table 1). 



2.2 SDSS cluster correlations 

We estimate the redshift space correlation function us- 
ing the estimator described by Hamilton (1993): 



&,(»•)= 4 
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where j = 1,2 and Ndd is the number of cluster pairs 
in the interval [r— Ar,r+Ar]. While, (Nrr) and (Ndr) 
is the average, over 10000 random simulations with the 
same properties as the real data (boundaries and red- 
shift selection function), cluster-random and random- 
random pairs, respectively. The random catalogues were 
constructed by randomly reshuffling the angular coordi- 
nates of the clusters (within the limits of the catalogue) , 
while keeping the same redshifts and thus exactly the 
same redshift selection function as the real data. 

Note that in order to take into account the possible 
systematic effects (eg. fraction of high-z clusters missed 
by the finding algorithm due to SDSS magnitude limit) 
in the different cluster subsamples we generate random 
catalogs, utilized the individual distance distribution of 
each subsample and not the overall SDSS cluster selec- 
tion function. We compute the errors on £s. (r) from 100 
bootstrap re-samplings of the data (Mo, Jing & Borner 
1992). 

We apply the correlation analysis to the Si and 
S2 subsamples evaluating £sj (r) in logarithmic inter- 
vals. In Fig. 2, we present the estimated two point red- 
shift correlation function (dots), divided according to 
richness class; strong clustering is evident. The dashed 
lines correspond to the best-fitting power law model 
£s-(f) = {fo/r) 1 , which is determined by the standard 
minimization procedure in which each correlation 
point is weighted by its error -1 . The fit has been per- 
formed taking into account bins with r > 5 h~ Mpc in 
order to avoid the signal from small, non-linear, scales 
while we have used no upper r cut-off (due to our a" 1 
weighting scheme, our results remain robust by varying 
the upper r limit within the 25 to 100 h -1 Mpc range). 
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Table 1. Results of the correlation function analysis for clusters of the two richness class (Si and S2 samples). Errors of the 
fitted parameters represent 2cr uncertainties. Finally, the r Q has units of h^ 1 Mpc. 



Sample No. of clusters r 7 r (7 = 1.8) 

51 200 20. 7±|;| 7 = l.et°A r ° = ^-^-H 

5 2 524 9-7li;a 7 = 2.0t ;£ r a = 9.8±\l 



In Fig. 3 we present the iso-Ax 2 contours (where 
^X 2 = X 2 - Xmin) in the 7 - r a plane. The Xmin is 
the absolute minimum value of the x 2 - The contours 
correspond to la (Ax 2 = 2.30) and 2a (Ax 2 = 6.17) 
uncertainties, respectively. In the insert of Fig. 3 we 
show the variation of A^ 2 around the best fit, once we 
marginalize with respect to the other parameter, while 
in Table 1, we list all the relevant information. For the 
Si cluster subsample (Abell R > richness) the best fit- 
ted clustering parameters are r a = 20.7^3 g/i _1 Mpc and 
7 = 1.6lg 4 which are in very good agreement with the 
values r a = 20.6 ± 1.5fr _1 Mpc and 1.5 ± 0.2 derived by 
Peacock & West (1992)*. Results for the S 2 subsample 
(APM richness) r D = 9.7±\%hT l Mpc and 7 = 2.0±° 7 5 
can be compared with those obtained by Dalton et al 
(1994); Bahcall & West (1992) and recently, from Plio- 
nis & Basilakos (2002), based on the APM cluster cata- 
log. They found a somewhat greater correlation length 
r ~ 12 - 13 h,- 1 Mpc. We can further estimate an 
upper limit of the correlation length using the expres- 
sion between the r and the mean cluster separation 
of Bachall & Burgett (1986), as modified by Bahcall & 
West (1992): r ,s 2 ^ 0.4ds 2 ~ 14.2 h' 1 Mpc (see also 
Dalton et al. 1994 and Croft et al 1997). 

In order to directly compare the correlation lengths 
of the two subsamples, we fixed the correlation function 
slope to its nominal value of 7 = 1.8 and we found 
r a = 19.8^3 1 and r = 9.81J'?. respectively (see last 
column of Table 1.). It is clear that the correlation 
length increases with cluster richness, as expected from 
the well-known richness dependence of the correlation 
strength. 

Finally, we have investigated the isotropy of the 
clustering signal for both subsamples, by examining the 
radial and tangential component of the SDSS correla- 
tion function £(<?, ir), with it the line-of-sight separation 
and a the perpendicular component of the cluster sepa- 
ration r (cf. Efstathiou et al. 1992). We have used bins 
of 20 h^ 1 Mpc width and in Fig. 4 we present the £ (a, ir) 
for both subsamples. It is evident that the £ (a, ir) con- 
tours are elongated along the line-of-sight direction, ir, 
up to ~ 40 h^ 1 Mpc. However, we suspect that this is 
not an indication of systematic effects related to line-of- 
sight projections but rather due to the extremely small 
width of the survey area (2.4 degrees in the declination 
direction which corresponds to ~ 25 h^ 1 Mpc at a red- 
shift of ~ 0.25), a fact that gives predominance to su- 
perclusters elongated along the line-of-sight with respect 
to those in the perpendicular direction (see Jing, Plionis 

* The robustness of our results to the fitting procedure was 
tested using different bins (spanning from 10 to 20) and we 
found very similar clustering results. 
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Figure 3. Iso-Ax 2 contours in the 7-ro parameter space for 
the Si (continuous line) and S2 (dashed line) samples. In the 
insert we show the variation of A\ 2 around the best fit once 
we marginalize with respect to the other parameter. 

& Valdarnini 1992 for the effects of superclusters elon- 
gated along the line-of-sight). Below we investiate our 
suspesion using the Hubble volume ACDM simulation 
(cf. Frenk et al 2000). 

2.3 Testing the robustness of the SDSS 
cluster correlations 

In order to test whether it is possible to recover the true 
underline cluster correlations from a survey with the ge- 
ometrical characteristics, selection function and richness 
of the SDSS, we have used the ACDM Hubble volume 
cluster catalogues (Colberg et al 2000). As an example, 
we present in Fig. 5 the underline Si Tike cluster corre- 
lation function, estimated from the whole volume (con- 
tinuous line) and the mean of 6 mock Si SDSS cluster 
samples (which contain around 200 clusters each). The 
mean clustering length of the SDSS mock samples is 
r ~ 18.5 ft -1 Mpc while that of the underline cluster 
population is r — 19.2 h^ 1 Mpc. It is evident that the 
SDSS survey is adequate to recover the underline clus- 
tering signal, albeit with a scatter of cr(£)/£ ~ 0.3 at 
separations, for example of, r ~ 15 h^ 1 Mpc. 

Furthermore, we address the issue of the observed 
anisotropies along the line of sight (see Fig. 4) by search- 
ing whether mock obervers show similar Si-like cluster- 
ing elongations along either their it or a directions. We 
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Figure 4. The radial and tangential anisotropy of the two 
point correlation function of the Si SDSS cluster subsample 
(while in the insert we show the results of the S2 subsam- 
ple). The transitions between different shadings correspond 
to fixed values of £(<r, it) = 1, 0.75, 0.5. The lines corresponds 
to the expected contours, using £(<r, n) = (ro/(<j 2 +tt 2 )) 7 . 
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Figure 6. The TZ = 6,2/6,1 for 50 mock Si like SDSS 
surveys for which 6,2 or 6,1 i s larger than 0.4. Note that 
1Z > 1 values indicate anisotropics along the 7r direction, 
while the arrow shows the observed SDSS TZ value. The insert 
plot shows the amplitude of the anisotropies, with respect 
to the central value (6,l)- Note that the continuous line 
divides anisotropies along the it direction (right of the line) 
and anisotropies along the <r direction (left of the line). 




Figure 5. Comaprison of the underline ACDM Si-like clus- 
ter correlation function (continuous line) with the mean of 6 
SDSS mock samples. 

quantify this anisotropy by the ratio, 1Z, of the £(cr, it) 
in the bins (ir,<r) = (0-20, 20-40) h' 1 Mpc (hereafter 
6,2) and (tt, a) = (20 - 40, - 20) h' 1 Mpc (hereafter 
6,1). 

We have selected 114 independent mock SDSS sur- 
veys (by spanning the z coordinate axis of the simula- 
tion) and we have found that in the majority of the cases 
(~ 60%) the value of 1Z is larger than unity, indicating 
a predominance of anisotropies along the 7r direction, 
while only ~40% of the cases it is less than one, indi- 
cating anisotropies along the a direction. 

We then investigated the amplitude of these S\ 
clustering anisotropies by deriving the distribution of 
the 1Z values for those observers that see relatively high 
correlation values; 6,2 > 0.4 and £2,1 > 0.4 (Fig. 6). In 
total we find 50 such observers out of which 37 (74%) 




Figure 7. As in Fig. 6 but for the S2 sample. 



show elongation along their line-of-sight. Therefore, it 
is evident, also due to the tail towards large 1Z values, 
that there are systematic anisotropies along the line- 
of-sight, which of course could be only due to the geo- 
metric characteristics of the mock cluster distribution. 
Furthermore in the insert of figure 6 we plot a scat- 
ter diagram between the normalized, by the value of 
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£(<7,7r) at the first bin (ie., < a, it < 20 h^ 1 Mpc), 
values of £1,2 and £2,1- It is apparent that the observed 
SDSS value (filled point) is roughly consistent with of 
the simulation derived values, although it appears to 
be an extremum. This indicates that the major part of 
the observed anisotropy is indeed due to the geometri- 
cal characteristics of the survey, as we have anticipated 
in the previous subsection, however we cannot exclude 
that some small contribution from intrinsic systemat- 
ics effects, like projection effects (cf. Sutherland 1988), 
could be present. 

A similar analysis of the S2 richness cluster cor- 
relations has shown that the percentage of mock ob- 
servers having significant clustering anisotropies along 
the 7r direction is significantly higher than those seeing 
anisotropies along the a direction (9 out of 10), however 
the amplitude of these anisotropies appear to be lower 
than in the observed case. In Fig. 7 we show the corre- 
sponding 1Z distribution and the scatter plot, in which 
it is evident that the amplitude of the true SDSS clus- 
tering anisotropy in the it direction is quite larger than 
what expected due to the survey geometrical character- 
istics. 

Therefore, we conclude that in the case of the Si 
cluster correlations there is no significant evidence for 
contamination by projection effects, while in the case of 
the S2 correlations we do have such indications. How- 
ever, in order to perform an a posteriori correction of 
£(r) (cf. Efstathiou 1992), for projection effects it would 
be necessary to disentangle first the effects of the sur- 
vey geometry, a task which at the present time is out of 
the scope of this work. Therefore, we caution the reader 
that all results based on the S2 cluster sample could be 
affected by the above mentioned systematic effect. 



3 MODEL CLUSTER CORRELATIONS 

It is well known (cf. Kaiser 1984; Benson et al. 2000) 
that assuming linear biasing the mass-tracer and dark- 
matter correlations, at some redshift 2, are related by: 

£modd(r, z) = £ D M(r, z)b 2 (z) , (2) 

where b(z) is the bias redshift evolution function. In the 
present work we have used the so called test particle 
bias model described by Nusser & Davis (1994), Fry 
(1996) and Tegmark & Peebles (1998). In this case the 
evolution of the correlation bias is developed assuming 
that only the test particle fluctuation field is related 
proportionally to that of the underling mass. Therefore, 
the bias factor as a function of redshift is 

*M = i + %^. (3) 

with 60 being the bias at the present time and D(z) 
the linear growth rate of clustering (described in sec- 
tion 4). It has been found (Bagla 1998) that, in the in- 
terval < 2 < 1, the above formula represents well the 
evolution of bias. Furthermore, the more accurate lin- 
ear bias evolution model given by Basilakos & Plionis 
(2001; 2003) is also very similar to the model of eq.(3) 
within 2 < 1. 



We quantify the evolution of clustering with epoch 
presenting the spatial correlation function of the mass 
£dm(t", 2) as the Fourier transform of the spatial power 
spectrum P(k): 

£ DM (r, 2) = D 2 (z)^Lj™ k 2 P(kf-^±dk , (4) 

where k is the comoving wavenumber. 

As for the power spectrum, we consider that of 
CDM models, where P{k) « k n T 2 {k) with scale- 
invariant (n = 1) primeval inflationary fluctuations. 
We utilize the transfer function parameterization as in 
Bardeen et al. (1986), with the approximate corrections 
given by Sugiyama's (1995) formula: 

TW = ln(1 2 +^ 49) [l + 3.89 9 + ( 16.1 g f + 
(5.46g) 3 + (6.7b?) 4 n 1/4 . 

with 



y n h 2 e xp [-n b - (2h) 1 / 2 n b /n \ y 1 

where k = 2tt/\ is the wavenumber in units of h Mpc -1 
and flj, is the baryon density. 

In the present analysis we consider flat models 
with cosmological parameters that fit the majority of 
observations, ie., fi m + Q,q = 1, H a = 100/ikm s _1 
Mpc -1 with h ~ 0.7 (cf. Freedman et al. 2001; Plionis 
2002; Peebles and Ratra 2002 and references therein), 
baryonic density parameter Qi,h 2 ~ 0.02 (e.g. Olive, 
Steigman & Walker 2000; Kirkman et al 2003) and a 
CDM shape parameter T = 0.17. In particular, we inves- 
tigate 3 spatially flat low-O m = 0.3 cosmological models 
with negative pressure and values of w = — 1 (ACDM), 
w = -2/3 (QCDM1) and w = -1/3 (QCDM2). Note 
that all the cosmological models are normalized to have 
fluctuation amplitude, in a sphere of 8 /i _1 Mpc radius, 
of cr 8 = 0.50(±0.1)fi m 7 (Wang & Steinhardt 1998) with 
7 = 0.21 -0.22w + 0.33O m . 



4 THE LINEAR GROWTH RATE OF 
CLUSTERING 

For homogeneous and isotropic flat cosmologies, 
driven by non relativistic matter and an exotic fluid 
(quintessence models) with equation of state, pq = wpQ 
and — 1 < w < 0, the Friedmann field equations can be 
written as: 

H 2 =(^) 2 = ^(p m + W ) (6) 
and 

- = -4ttG w+- )PQ + -Pm] , (7) 
a 3 3 

where a(t) is the scale factor, p m oc a(t)~ 3 is the mat- 
ter density and pQ oc a(t)~ 3( - 1+w ^ is the dark energy 
density. 

The time evolution equation for the mass density 
contrast, modeled as a pressureless fluid has general so- 
lution of the growing mode (Peebles 1993): 
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D + 2H(t)t) = 47rGp m D 



(8) 



where dots denote derivatives with respect to time. 
From the equations describing the Friedmann model, 
it follows that H + H 2 = -4ttG[(w + l/3)p Q + (l/3)p m ]. 
Differentiating this relation and using p m = — 3ffp m 
p'q = —3(1 + w)HpQ we obtain 

H + 2HH = 4ttG(1 + w)(w + \)p Q H + 4nGp m H . (9) 

o 

Therefore, it turns out that if w = — 1 (ACDM) ortu = 
-1/3 (QCDM2) then H(t) is a decaying mode of eq.(8). 
In that case, the growing solution (Peebles 1993) as a 
function of redshift is: 



D{z) = 



(1 + x) 



da; 



2 7 z 

where we have used the following expressions: 

e(z) = [n m (i + «) 3 + n (i + ^) 3(1+t " ) ] 1/2 

dt 1 

dz ~ ~~ H E(z){\ + z) ' 



(10) 

(11) 
(12) 



The Hubble parameter is given by: H(z) = H E(z), 
while f2 m = 8ivGpo/3Ho (density parameter), Qq = 
8ttGpq/3H* (dark energy parameter) at the present 
time, which satisfy fl m + Qq = 1 and finally H is the 
Hubble constant. In addition to to m {z) also Qq(z) could 
evolve with redshift as 

On(l + z) 3 



Q m (z) = 

and 



E 2 (z) 



fig(i + *) 



3(1+™) 



(13) 



(14) 



It is interesting to mention that in a flat low-O m 
with w = —1/3 model, the equation of state pq = 
— (1/3) pg leads to the same growing mode as in an open 
universe, despite the fact that this quintessence model 
has a spatially flat geometry! Therefore, as the time 
evolves with redshift, utilizing equations (12), (11) and 
the relation 



o tt2 

47rGp m = ^a m (l + «) 3 



(15) 



then the basic differential equation for the evolution of 
the linear growing mode takes the following form: 



with basic factors, 

1 1 AE(z) 



P(z) 
and 

Q(z) 



l + z 



+ 



E(z) dz 
2E 2 (z) 



(16) 



(17) 



(18) 



We find that eq.(16) has a decaying solution of the form 
Di(z) = (1 + z) n only for w = -2/3, with n = 3/2. 
The second independent solution of eq.(16) can be found 
easily from the following expression: 
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Figure 8. The variance A\ 2 around the best fit b value 
for various cosmological models. The left and the right panel 
correspond to the S2 and Si samples, respectively. 



f°° dx r r 

D(z) = D 1 (z) j ^) exp - j P(t)dt 
s to the folio 

. r ± 

h (1+-) 



(19) 



which finally leads to the following growing mode: 
D(z) = (l + zf' 2 I 



) 2 E{x) 



5 THE SDSS CLUSTER BIASING 

In order to quantify the cluster bias at the present time 
we perform a standard x 2 minimization procedure (de- 
scribed before) between the measured correlation func- 
tion of the SDSS galaxy clusters with those expected in 
our spatially flat cosmological models 

£ Si M-£nodelMo) n2 



X{bo) 



E 



(21) 



where a % is the observed correlation function (boot- 
strap) uncertainty. 

In Fig. 5 we present, for various cosmological mod- 
els, the variation of A\ 2 = X 2 (^°) — Xmin(^o) around the 
best 60 fit, for the different richness class (left panel for 
S2 and right panel for Si). 

To this end, owing to the fact that the observational 
data are analyzed in redshift space, the correlations 
should be amplified by the factor K(/3) = 1+2/3/3+/3 2 /5 
(Hamilton 1992) where /3 ~ fC/bo. We utilize the 
generic expression for a, defined by the Wang & Stein- 
hardt (1998): 

! ^ ^'".'-y d-ft.) -(22) 



5 - w/(l - w) 125 (1 - 6w/5) 3 

In Table 2 we list the results of the fits for our two cluster 
catalogs, ie., the cosmological models and the value of 
the cluster optical bias, b , at the present time, as well 
as the redshift distortion /3 parameter and a measure 
of the K((3) correction. We find that the redshift space 
distortions effect increases £sj (r) by a factor of ~ 12 — 
26%. 

In Fig. 6, we plot the measured £sj (r) (filled sym- 
bols) of our two samples with the estimated two point 
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Table 2. List of the theoretical clustering model fitting analysis. Note that the errors of the fitted parameters represent 2a 
uncertainties. 



Index b /3 K((3) 



ACDM-Si 2A+° 2 3 0.20+°;°^ 1-14 j^'.! 

QCDMl-Si 2.6±°;^ 0.19±g;g|^ 1.13±%' t 

QCDM2-Si 2.7l°;^ 0.18±g;gf| 1.121°-; 

ACDM-S2 l-4loi 0.35+°'°23 l-26±g-; 

QCDMI-S2 l-7±° i 0.29j;° : °^ 1.21±g : ; 
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Figure 9. Comparison of the observed and model SDSS clus- 
ter correlation function: Left panel, S2 sample; right panel, 
Si sample. The observational data are represented by filled 
symbols. 



correlation function for all three cosmological models. 
We should conclude that the behavior of the observed 
two point correlation function of the galaxy clusters is 
sensitive to the different cosmologies with a strong de- 
pendence on the present time bias. By separating be- 
tween low and high richness regimes, we obtain results 
being consistent with the hierarchical clustering sce- 
nario, in which the rich clusters are more biased tracers 
of the underlying matter distribution with respect to 
the low richness clusters. 

We can put some further cosmological constraints, 
comparing our clustering results with those based on 
large-scale dynamics. For example Branchini & Plionis 
(1996) using the cluster dipole after reconstructing the 
spatial distribution of Abell/ACO R > clusters found 
/?Abeii = 0.21 ± 0.03. Also, Branchini at al. (2000) com- 
paring the density and velocity fields of the Abell/ACO 
cluster distribution with the corresponding POTENT 
fields (using the MARK III galaxy velocity sample) , ob- 
tained /^potent = 0.22 ± 0.08. Comparing the latter (3- 
results with our clustering predictions (Table 2) we can 
conclude that for the Si sample (Abell R > richness) 



the only model which fails (although marginally) to re- 
produce the large-scale dynamical results is QCDM2 
(w = -1/3). 



6 CONCLUSIONS 

We have studied the clustering properties of the SDSS 
galaxy clusters in redshift space. We have divided the 
total sample in two richness subsamples; roughly corre- 
sponding to Abell R > (A^gai > 30 members) and to 
APM (N gai > 20 members) clusters. We find that if the 
two point cluster correlation function is modeled as a 
power law, £(r) = (r /r) 7 , then the best-fitting parame- 
ters are (a) r D = 20.7±|;| h" 1 Mpc with 7 = 1.6±^t and 
(b) r D = 9.7±i J h' 1 Mpc with 7 = 2.0±° l respectively. 
We have also found that the Abell-like sample is not sig- 
nificantly affected by projection effects, and its appar- 
ent clustering elongation along the line-of-sight is due 
to the survey geometry. However, the APM-like sample 
appears to be somewhat affected by projection effects, 
showing a clustering elongation along the line-of-sight 
larger than what expected from the survey geometry. 

Comparing the cluster correlation function with the 
predictions of 3 spatially flat quintessence models (hav- 
ing Q m = 0.3), we estimate the cluster redshift space 
distortion parameter K(/3) and we conclude that the 
amplitude of the cluster redshift correlation function 
increases by a factor of ~ 12 — 26% (depending on the 
richness class) . Finally, comparing our clustering results 
with those of dynamical analysis, based on the large 
scale motions, we find that the flat cosmological mod- 
els with w < —0.6 are consistent with the observational 
results. 
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